Part of the D - dimensional Spiked harmonic 

oscillator spectra 



Omar Mustafa and Maen Odeh 
Department of Physics, Eastern Mediterranean University 
G. Magusa, North Cyprus, Mersin 10 - Turkey 
email: omustafa.as@mozart.emu.edu.tr 

Abstract 

The pseudoperturbative shifted - I expansion technique PSLET 
[5,20] is generalized for states with arbitrary number of nodal ze- 
ros. Interdimensional degeneracies, emerging from the isomorphism 
between angular momentum and dimensionality of the central force 
Schrodinger equation, are used to construct part of the D - dimen- 
sional spiked harmonic oscillator bound - state spectra. PSLET re- 
sults are found to compare excellently with those from direct numerical 
integration and generalized variational methods [1,2]. 
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1 Introduction 



The simplest model of realistic interaction potentials in atomic, molecular, 
and nuclear physics is provided by the spiked harmonic oscillator 

V(q) = c iq 2 + c 2 q- b , ci,C2,6>0 , ge(0,oo). (1) 

The construction of its bound - states has attracted attention over the last 
few years [1-9]. It is an interesting model not only because of being a singu- 
lar potential representing a repulsive core in realistic interactions, but also 
because of its intrinsic properties in view of mathematical physics [10-16]. 
However, most of the studies on this model potential (1) were devoted to 
one spatial dimension ( ID, the hyperquantum limit in view of Herschbach 
[17,18]). It was just very recently, to the best of our knowledge, that Hall and 
Saad have generalized their variational analysis [la] and smooth transforma- 
tion [2] methods, VAM and STM, respectively, to the D - dimensional case 
and studied its bound - states. They have also used direct numerical integra- 
tion (DNI) for comparison purposes. It is therefore interesting to carry out 
systematic studies of the bound - state spectra generated by this interesting 
class of singular potentials (1). 

On the other hand, results from exactly solvable potentials ( an inter- 
esting field of mathematical physics in itself) are essential ingredients for 
the description of realistic physical problems [1-5,19]. The solutions of these 
can be used in perturbation and pseudoperturbation theories, or they can 
be combined with numerical calculations. Nevertheless, in the simplest case, 
analytical calculations can aid numerical studies in areas where numerical 
techniques might not be safely controlled. For example, when bound - state 
wave functions with arbitrary nodal zeros are required for certain singular 
potentials (a next level of complexity), analytical solutions can supply a 
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basis for numerical calculations. Moreover, in many problems the Hamil- 
tonian does not contain any physical parameter suitable for a perturbation 
expansion treatment. More often, the Hamiltonian contains physical param- 
eters, but, typically, zeroth - order solutions for special values of these are 
not tractable or good starting approximations. One would therefore resort 
to variational calculations [1], pseudoperturbation expansions ( artificial in 
nature) [5,18-26], etc. 

Recently, we have introduced a pseudoperturbative shifted - I (lis the an- 
gular momentum quantum number) expansion technique ( PSLET) to solve 
for nodeless states of Schrodinger equation. It simply consists of using 1/7 
as a pseudoperturbation parameter, where I — I — (5 and (3 is a suitable shift. 
The shift (3 is vital for it removes the poles that would emerge, at lowest 
orbital states with /=0, in our proposed expansions below. Our analytical, 
or often semianalytical, methodical proposal PSLET has been successfully 
applied to quasi - relativistic harmonic oscillator [20], spiked harmonic oscil- 
lator [5], anharmonic oscillators [21], and to the two - dimensional (flatland, 
in view of Godson and Lopez - Cabrera in [17]) hydrogenic atom in an arbi- 
trary magnetic field [22]. 

Encouraged by its satisfactory performance in handling nodeless states, 
we generalize PSLET recipe ( in section 2) for states with arbitrary number 
of nodal zeros, k > 0. Moreover, in the underlying "radical" time - indepen- 
dent radial Schrodinger equation, in h — m — 1 units, 



the isomorphism between orbital angular momentum I and dimensionality D 
invites interdimensional degeneracies to obtain [17]. Which, in effect, allows 
us to generate the ladder of excited states for any given k and nonzero I from 
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the 1=0 result, with that k, by the transcription D — > D + 21. That is, if 
Ekj(D) is the eigenvalue in D - dimensions, then 

E H {2) = £ M _i(4) = ■■■ = E ktl (2l) = E kfi (2l + 2) (3) 
for even D, and 

E ktl (3) = £ M _i(5) = • • • = E ktl (2l + 1) = E kfi (2l + 3) (4) 

for odd D. For more details the reader may refer to ref.s [17,18,27]. We 
therefore calculate, in section 3, the energies for D = 2 and D = 3 spiked 
harmonic oscillators, for a given number of nodes k and different values of I, 
and construct part of its D - dimensional bound - state spectra. We compare 
our results with those reported by Hall and Saad via generalized variational 
analysis VAM, and direct numerical integration DNI methods [1,2]. Section 
4 is devoted for concluding remarks. 

2 The generalization of PSLET 

With the shifted angular momentum, equation (2) reads 

where Q is a constant that scales the potential V(q) at large - Id limit ( the 
pseudoclassical limit [17]) and is set, for any specific choice of Id and k, equal 
to I 2 at the end of the calculations. Here l D — I + (D — 3)/2, to incorporate 
the interdimensional degeneracies associated with the isomorphism between 
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angular momentum and dimensionality D. Hence, / — > I — Id — ft through 
out this paper. Next, we shift the origin of the coordinate system through 
x = l l l 2 (q — q )/Qo, where q Q is currently an arbitrary point to be determined 
below. Expansions about this point (see Appendix for more details), x = 
(i.e. q = q Q ), obviously localize the problem at an arbitrary point q Q and 
the derivatives, in effect, contain information not only at q Q but also at any 
point on g-axis, in accordance with Taylor's theorem. It is then convenient 
to expand E k ; as 



E k , = E eS r». 



(6) 



Equation (5) thus becomes 



1 rl 2 00 

2 dx 2 ^ 



n=l 



**,!(*)■ (7) 



Up to this point, one would conclude that the above procedure is nothing 
but an imitation of the eminent shifted large-N expansion (SLNT) [25,26,28- 
30]. However, because of the limited capability of SLNT in handling large- 
order corrections via the standard Rayleigh-Schrodinger perturbation theory, 
only low-order corrections have been reported, sacrificing in effect its precise- 
ness. Therefore, one should seek for an alternative and proceed by setting 
the wave functions with any number of nodes as 

*k,i( x (q)) = F kA x ) exp{U Kl {x)). (8) 



In turn, equation (7) readily transforms into the following Riccati equation: 
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F k ,i(x) 



1 oo 

K» + u' Kl {x)u' Kl {x)) + £ ^ (n) (*)r n/2 



n=0 



n=l 



(9) 



where the primes denote derivatives with respect to x. It is evident that this 
equation admits solution of the form 



00 (JXJ 

KM = E ui n \x) + E 4 n) (*) r^ 2 , (10) 



?1=0 



71=0 



OO fc— 1 



n=0p=0 



(11) 



where 



rt+l 



/ j Dm,n,k 

m=0 



2m- 1 



; £W = 0, (12) 



n+l 



X 



2m 



m=0 



Substituting equations (10) - (13) into equation (9) implies 



(13) 



F k ,i(x) 



~2 E (^ n) ^~ n/2 + G i n) r {n+1)/2 J 



1 oo n 

c\ / j / j \ h h k k 



n=0 m=0 
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OO CO 

+ 2ul m) G { r rn) i~^/ 2 ) + ]T v^r n ' 2 - J2 qlE^r™ 



n=0 n=l 
rWf-n/2 , ri(»)T-(i»+l)/2 



^^(n) / -- ri/2 + G (n) [ _(„ + l)/2^ 



.n=0 



- ^» = (14) 



The above procedure obviously reduces to the one described by Mustafa and 
Odeh [5,20-22], for k = 0. Moreover, the solution of equation (14) follows 
from the uniqueness of power series representation. Therefore, for a given k 
we equate the coefficients of the same powers of I and x, respectively. For 
example, when k = 1 one obtains 



01,0,1 = -w, U{°\x) = - w x, 



(15) 



Ci, ,i - 



a (1) - 
"o,i — 



Co,o,: 



w 



(16) 



C 0) o 1 i = -(2C 1)0l i + 2/3 + 1), 



(17) 



02,2,1 = ^ - 5 4 ) , 



(18) 



0i,2,i = ^ (J 02, 2 ,i + Cb.0,1 C 1>0>1 - \{2(5 + 1) ) , 



(19) 



E (o) = 1 W + l) (i) c _ 3 0i, 2 ,i _ CV 
i,z n 2 1 2 ^ -'-'0'-'- 2 2 



9o 



(20) 
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etc. Here, we reported the nonzero coefficients only and give the definitions 
of the related parameters in the Appendix. One can then calculate the en- 
ergy eigenvalues and eigenfunctions from the knowledge of Cm,^, Dm,n,k, 
and a^l in a hierarchical manner. Nevertheless, the procedure just described 
is suitable for a software package such as MAPLE to determine the energy 
eigenvalue and eigenfunction corrections up to any order of the pseudoper- 
turbation series (6). 

Although the energy series, equation (6), could appear divergent, or, at 
best, asymptotic for small I, one can still calculate the eigenenergies to a very 
good accuracy by forming the sophisticated [N,M] Pade approximation [24] 

P#(l/J) = (P + PJI + ■■■ + Pm/T M )/{1 + Qi/l + ■■■ + q N /l N ) 
to the energy series (6). The energy series (6) is calculated up to EfJ /I 8 by 

E k ,i = PE™ + + ■■■ + E$/T> + (21) 
and with the P£(l/l) Pade approximant it becomes 

E k MA} = PEl; l 2) + p*(i/i). (22) 

Our recipe is therefore well prescribed. 

3 D - spiked harmonic oscillator spectra 

In this section we consider the spiked harmonic oscillator potential (1) and 
illustrate the above mentioned procedure. The substitution of equation (1) 
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in (45), for k > 0, implies 



w 



8c 1 q + bc 2 (b-2)q~ (b+1) 
\ 2ciq Q - bc 2 qo ih+1) 



, /3 = --(l + [2h + l]w). (23) 



Equation (44), in turn, reads 



l D + - I 1 + [2k + 1] 



8cig D + bc 2 (b - 2)q 
\ 2c x q Q - bc 2 qo {b+1) 



(6+1) 



2 4 L _ ^ -(&+2) 

y y c i 2 g ° ' 



(24) 



which is explicit in g D . However, in the absence of a closed - form solution 
for q Q , which is often the case ( hence the notion that PSLET is often semi- 
analytical), numerical solutions of (24) could resolve this issue. Once q is 
determined the coefficients C mtTlj k, D m>n) k, and a^l are determined in a se- 
quential manner. Hence, the eigenvalues, equation (21), and eigenfunctions, 
equations (10)-(13), are calculated in the same batch for each value of k, D, 
I, Ci, c 2 , and b. 

Table 1 shows PSLET results for the ground - state energies, covering a 
wide range of the coupling c 2 when b = 2.5, along with those reported by 
Hall and Saad [la], via a generalized variational analysis and direct numerical 
integration methods. Using the interdimensional degeneracies, equations (3) 
and (4), or directly the dimensionality D in l D , we display the energies for 
V(q) = (<? 2 + 10/g 19 )/2 in table 2. Clearly, our results compare excellently 
with those from direct numerical integrations. However, it should be noted 
that in [5] we have calculated the energy series up to E^J /I 4 correction. 
Therefore, slight discrepancies obtain between the present results in table 1 
and those reported in table 2 of [5]. 

Adhering to the implicated wisdom in equations (3) and (4), that the 
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two - and three - dimensional ( 2D and 3D, respectively) cases are the basic 
ingredients of the energy ladder at larger dimensions, we report (in table 3) 
the 2D - and 3D - nodal bound - state energies when the coupling c 2 = 1000 
and b = 0.5, 1, • • • , 2.5, 3. The stability of the last three approximants of the 
Pade sequence indicate that the results are exact. For more details on this 
issue the reader may refer to ref.s [20,24]. Nevertheless, our results E 0:0 for 
the 3D - spiked harmonic oscillator are in exact accord with those from direct 
numerical integrations [2] . Following the same strategy, we display in table 4 
the k — 1 and 2 nodal bound - state energies for V(q) = (q 2 + 1000/g 3/2 )/2. 
Eventually, the leading term of PSLET, PeI~ 2 \ turns out to be a good 
starting approximation. Tables 1,2 and 4 bear this out. 

Moreover, for the spiked harmonic oscillator, with 6 = 2, one would 
rewrite the effective potential term (1(1+ 1) + c 2 )/2g 2 + g 2 /2 as Z'(/' + l)/2g 2 + 
q 2 /2 with I' = -1/2 + + 1/2) 2 + c 2 . For this particular case, PSLET pro- 
cedure yields, respectively, w = 2, (3 = -(2k + 3/2), 1 = 2k + l' + 3/2, q 2 Q = I, 
t 1 E <y ~^ > = 2k + l' + 3/2 ( the exact well known energies), 

<? = ^; = - = ^ = - = ^=0, (25) 
and when k — 0, for example, 



TJ , x i ( y 2 . y 3 y\ y 5 y 6 . y 7 y 8 . \ 
U ^ x) = -2{ y -Y + J-J + J-J + T-J + ---) 

.r( y 2 . y 3 y A . y 5 y 6 . y 7 y s . \ 
y - h, (26) 



where y = xl 1//2 . Obviously, the terms in brackets in equation (26) are the 
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infinite geometric series expansions for ln(l+y). Equation (26) thus becomes 



U o/ (x) = ln(l + y)- 1 ' 2 + ln(l + y) 1 - ly - 



(27) 



Hence Eq.(8) ( with F ol >(x) = 1 from (11)) reads 



/+1 p -9 2 /2 



(28) 



the exact well known solutions [31], where N ol > are the normalization con- 
stants. Proceeding exactly as above, one could obtain the well known so- 
lutions with k > 1. However, this already lies far beyond the scope of our 
present methodical proposal. 

Hall and Saad [la] have therefore used, indirectly, the transformation of 
the angular momentum quantum number and cast the Hamiltonian of the 
spiked harmonic oscillator (1) as 



Where Ih — — 1/2 + y (I + 1/2) 2 + A, and A is used as a further variational 
refinement in their generalized variational analysis method. They found that 
A = C2 is a good general estimate for the value of A. Indeed this optimum 
value of A, which reduced substantially the number of the basis function 
needed for a given accuracy in [la], enhances the convergence and accuracy 
of approximation methodical recipes. Practically, it minimizes the effect 
of the perturbation term C2q~ b over the harmonic oscillator one (with the 
irrational quantum number /'), especially for values of b — ► 2. In table 5, 
the results of PSLET are obtained using such prescription. They compare 




(29) 
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excellently with direct numerical integrations and do not contradict with the 
upper bounds from the generalized variational estimates. 

4 Concluding remarks 

We have generalized our pseudoperturbative shifted - I expansion technique 
PSLET [5,20-22] for states with arbitrary number of nodal zeros, k > 0. 
Starting with the central force problem, represented by the radial Schrodinger 
equation, and augmenting the orbital angular momentum by I — > Id = 
I + (D — 3)/2, we have incorporated interdimensional degeneracies. To test 
PSLET performance, we have treated the spiked harmonic oscillator prob- 
lem in D - dimensions, and used results from direct numerical integrations 
and generalized variational analysis methods [1,2] to compare with. The 
comparison is readily satisfactory. 

The salient features of the attendant proposal PSLET are in order. 

It avoids troublesome questions as those pertaining to the nature of small 
parameter expansions, the trend of convergence to the exact numerical val- 
ues ( marked in tables 1-3 and 5), the utility in calculating the eigenvalues 
and eigenfunctions in one batch to sufficiently higher - orders ( documented 
through the solution (28) of (1), with 6 = 2), and the applicability to a wide 
range of potentials. Provided that the potential V(q) gives rise to one mini- 

(—2) 

mum of E kl and an infinite number of bound - states. Moreover, beyond its 
promise as being quite handy ( on the computational and practical methodi- 
cal sides), it offers a useful perturbation prescription where the zeroth - order 
approximation l 2 E\r^ inherits a substantial amount of the total energy. 

The above has been a very limited review and a number of other useful 
approaches such as those presented by Papp [9,32], Bender and Wu [33], etc, 
have been left unattended. However, their accomplishments are indeed of 
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actual novelties. 

Finally, the scope of PSLET applicability extends beyond the present D - 
dimensional spiked harmonic oscillator model. It could be applied to angular 
momentum states of multi - electron atoms [34-36], relativistic and non - 
relativistic quark - antiquark models [37], etc. 
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Appendix 



Although some of the following expressions have appeared in previous 
articles [5,20-22], we would like to repeat them to make this article self con- 
tained. 

Expansions about x = (i.e.q = q Q ), yield 



1 



= E(-iy 



n=0 



Q 2 o 



X ' 



(30) 



v{x( q )) = 1 i d ^A ^ r^. 



n=0 



(31) 



Equation (5) thus becomes 



Id 2 a 2 ~ 



(32) 



with 



I 



V(x(q)) = q 2 l 



1 V(go) 



+ (2/3 + !)£(-!)" 



n=l 

oo 



n=0 



.(" + !), 
2 

>+l) 



U n r n/2 + £ B n x n z" 



(n-2)/2 



n=3 



-;r 



n^-(n+2)/2 



(33) 



-Bti 



-l) n 



1) 



+ 



v d(Q J n\Q 



(34) 
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Equation (32), along with (33) and (34), is evidently the one - dimensional 
Schrodinger equation for a perturbed harmonic oscillator 



X k (x) = \ k X k (x), 



(35) 



where w 2 = 2B 2 , 



£ = I 



1 J% 

2 Q 



2/3+1 /3Q3 + 1) 
+ 2 + 21 



(36) 



and P(x) represents the remaining terms in eq.(33) as infinite power series 
perturbations to the harmonic oscillator. One would then imply that 



+ 



1 Q 2 qV(Qq) 

2 Q 

W+l) + A (o) 



2/3 + 1 1. 

oo 

+ e A^r-«, 

n=2 



(37) 



and 



n=-2 



(38) 



Hence, equations (37) and (38) yield 



~2^ + Q 



(39) 



^k,i 



'2(3 + 1 1. 



(40) 
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E\ 



(0) _ i rw+i) . x(o) 



(41) 



'fc,Z — 2 O ^ ^ 



o 



4? = Ai n) /? 2 ; 



n > 1. 



(42) 



Where q Q is chosen to minimize E). 



(-2) 



dE, 



.(-2) 



= and 




> 0. 



(43) 



dg 



Hereby, V(q) is assumed to be well behaved so that E { k l has a minimum g D 
and there are well - defined bound - states. Equation (43) in turn gives, with 



Consequently, the second term in Eq.(33) vanishes and the first term adds a 
constant to the energy eigenvalues. It should be noted that the energy term 



momentum L z =l executing circular motion of radius q Q in the potential V(q ). 
It thus identifies the zeroth - order approximation, to all eigenvalues, as a 
classical approximation and the higher - order corrections as quantum fluc- 
tuations around the minimum q , organized in inverse powers of I. The next 
correction to the energy series, lE^T^, consists of a constant term and the 
exact eigenvalues of the harmonic oscillator w 2 x 2 /2.The shifting parameter 
(3 is determined by choosing lEj~^^—0. This choice is physically motivated. 




(44) 




corresponds roughly to the energy of a classical particle with angular 
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In addition to its vital role in removing the singularity at I — 0, it also re- 
quires the agreements between PSLET eigenvalues and eigenfunctions with 
the exact well known ones for the harmonic oscillator and Coulomb poten- 
tials. Hence 



P = - 



~ 2 + (k + -)w 



(45) 



where w = y3 + q V" (q )/V (q Q ), and primes of V(q ) denote derivatives 
with respect to q Q . Then equation (33) reduces to 



I 



V(x(q)) = q£l 



1 V(Qq) 
2q 2 o Q 



+ $> (n) (z)r n/2 , 



n=0 



(46) 



where 



(47) 



v W(x) = -{2/3 + l)x + B 3 x 3 , (48) 

and for n > 2 



vW(x) = B n+2 x n+2 + (-l) n (2/3+1) 



(!L+1) x n 



+ -1 



- 1) x(™- 2 ). 



(49) 
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Table 1: 3D ground - state energies, in h — m — 1 units, for V(q) = 
(q 2 + c 2 /q 5/2 )/2. Where E P represents PSLET results, Eq.(41), and l 2 E^ 
is its zeroth - order approximation. E[4, 4] shows the effect of the P±(l/l) 
Pade approximant ,Eq.(42). E V am from VAM, and E DNI from DNI [la]. 



c 2 


p£(-2) 


Ep 


£[4,4] 


EyAM 


Edni 


1000 


44.003142 


44.9554848 


44.9554848 


44.955485 


44.955485 


100 


16.666664 


17.541890 


17.541890 


17.541890 


17.541890 


10 


7.00149 


7.73515 


7.73510 


7.73511 


7.73511 


1 


3.84771 


4.31578 


4.31413 


4.32326 


4.31731 


0.1 


3.11132 


3.26984 


3.26633 


3.29602 


3.26687 


0.01 


3.0116 


3.0341 


3.0344 


3.0392 


3.0367 


0.001 


3.0012 


3.0035 


3.0040 


3.0041 


3.0040 
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Table 2: D — 2, • • • , 10 ground - state energies, in h — m — 1 units, for 
the potential V(q) = (q 2 + 10/g L9 )/2. Where E P represents PSLET results, 
Eq.(41), and PE 1 -^ is its zeroth - order approximation, #[4,4] shows the 
effect of the P%{l/T) Pade approximant , Eq.(42). E VAM from VAM, and 
E DNI from DNI [la]. 



D 




Ep 


#[4,4] 


EyAM 


Edni 


2 


7.581 139 


8.485 461 


8.485 369 


8.485 384 


8.485 378 


3 


7.919 880 


8.564 352 


8.564 355 


8.564 358 


8.564 356 


4 


8.339 920 


8.795 436 


8.795 440 


8.795 440 


8.795 440 


5 


8.840 678 


9.163 092 


9.163 093 


9.163 093 


9.163 093 


6 


9.416 352 


9.646 701 


9.646 701 


9.646 701 


9.646 701 


7 


10.058 042 


10.225 045 


10.225 045 


10.225 045 


10.225 045 


8 


10.755 870 


10.879 077 


10.879 077 


10.879 077 


10.879 077 


9 


11.500 402 


11.592 982 


11.592 982 


11.592 982 


11.592 982 


10 


12.283 349 


12.354 183 


12.354 183 


12.354 183 


12.354 183 
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Table 3: 2D - and 3-D - nodeless states energies, with / = 0,---,4 (in 
h — m — 1 units), for the potential V(q) = (q 2 + 1000/V)/2. Where E j 
represents PSLET results with the P%{1/1) Pade approximant, Eq.(42). 



D 


b 


-E'0,0 


Eo,i 


E>0,2 


^0,3 


Eo,4 


2 


0.5 


415.886751 


415.898889 


415.935293 


415.995938 


416.080780 




1 


190.719321 


190.735267 


190.783089 


190.862739 


190.974135 




1.5 


104.404517 


104.427341 


104.495769 


104.609681 


104.768874 




2 


65.245553 


65.277168 


65.371918 


65.529521 


65.749510 




2.5 


44.945030 


44.986838 


45.112071 


45.320150 


45.610129 




3 


33.303511 


33.356491 


33.515080 


33.778229 


34.144222 


3 


0.5 


415.889786 


415.914059 


415.962588 


416.035338 


416.132258 




1 


190.72331 


190.755196 


190.818940 


190.914475 


191.041704 




1.5 


104.41022 


104.455860 


104.547051 


104.683633 


104.865367 




2 


65.253459 


65.316665 


65.442888 


65.631753 


65.882705 




2.5 


44.95549 


45.039054 


45.205805 


45.454976 


45.785438 




3 


33.31676 


33.422634 


33.633677 


33.948503 


34.365078 
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Table 4: 2D and 3D k - state energies, in h — m — 1 units, for the potential 
V(q) = (q 2 + 1000/g 3 / 2 )/2. Where E P represents PSLET results, Eq.(41), 
and 1 2 E(~ 2 ^ is its zeroth - order approximation, -E[4,4] shows the effect of 
the Pf{l/l) Pade approximant, Eq.(42). 



D 


k 


I 




Ep 


£[M] 


2 


1 





105.40419 


108.15083 


108.15083 






1 


105.67466 


108.17379 


108.17379 






2 


105.96940 


108.24263 


108.24263 






3 


106.28970 


108.35721 


108.35721 


3 







105.53648 


108.15657 


108.15657 






f 


105.81892 


108.20248 


108.20248 






2 


106.12628 


108.29421 


108.29421 






3 


106.45983 


108.43160 


108.43160 


2 


2 





107.3876 


111.9017 


111.9017 






f 


107.7382 


111.9248 


111.9248 






2 


108.1127 


111.9940 


111.9940 


3 







107.5600 


111.9075 


111.9074 






f 


107.9224 


111.9536 


111.9536 






2 


108.3092 


112.0459 


112.0459 
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Table 5: k = 2 and 1 — 1,2 energies, in h — m — 1 units, for the potential 
V(q) = (q 2 + 10/g 2 - 1 )/2. Where E 2 ^ P represents PSLET results, Eq.(41), 
E 2 ,i,v from VAM, and E 2 ,i, ex from DNI [la]. E 2 ,/[4,4] shows the effect of the 
P£(l/l) Pade approximant, Eq.(42). 



D 


E 2 ,l,ex 


E 2 xv 


E 2j i : p 


E 2 ,i [4, 4] 


E 2)2> p 


E 2 , 2 [4, 4] 


2 


16.543629 


16.543648 


16.541951 


16.543627 


17.380817 


17.381708 


3 


16.904445 


16.904446 


16.903172 


16.904444 


17.954856 


17.955444 


4 


17.381708 


17.381709 


17.380817 


17.381708 


18.606695 


18.607067 


5 


17.955444 


17.955446 


17.954856 


17.955444 


19.320461 


19.320691 


6 


18.607067 


18.607070 


18.606695 


18.607067 


20.083266 


20.083406 


7 


19.320691 


19.320693 


19.320461 


19.320691 


20.884936 


20.885021 


8 


20.083406 


20.083407 


20.083266 


20.083406 


21.717556 


21.717608 


9 


20.885021 


20.885022 


20.884936 


20.885021 


22.574996 


22.575027 


10 


21.717608 


21.717608 


21.717556 


21.717608 


23.452505 


23.452524 
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